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Current LHC searches for new colored particles generally focus on their pair production channels 
and assume any single production to be negligible. We argue that such an assumption may be 
unnecessary in some cases. Inclusion of model dependent single productions in pair production 
searches (or vice versa) can give us new information about model parameters or better exclusion 
limits. Considering the example of the recent CMS search for first generation scalar leptoquarks in 
the pair production channel, we illustrate how single productions can be systematically included in 
the signal estimations and demonstrate how it can affect the mass exclusion limits and give new 
bounds on leptoquark-lepton-quark couplings. We also estimate the effect of the pair production in 
the more recent CMS search for scalar leptoquarks in single production channels. 
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I. INTRODUCTION 

Various beyond the Standard Model (BSM) theories 
predict a host of new heavy particles at the TeV range. 
The Large Hadron Collider (LHC) at CERN is search¬ 
ing for many of these particles. In a hadron collider like 
the LHC, a heavy colored particle could be produced via 
different mechanisms. Among them, usually the strong 
interaction mediated pair production dominates. Other 
direct production processes like single productions gen¬ 
erally depend on some model specific couplings (Ax). 
Hence, many present direct searches for new colored par¬ 
ticles at the LHC (and related phenomenological studies) 
focus on their pair productions and ignore single pro¬ 
ductions assuming the couplings controlling them to be 
smaller than the strong coupling. With this assumption 
the results (or predictions) seemingly become model in¬ 
dependent. But this, however, raises a question - is it 
justified to ignore single production completely? Obvi¬ 
ously, the question is relevant if any such ignored Ax is 
indeed not very small in the nature. But, even if Ax 
is small, cross-sections of some single productions could 
still be comparable to the pair production cross-section 
in some region of parameter-space. This is possible when 
the particle being produced is quite heavy so that the ex¬ 
tra phase-space suppression received by the pair produc¬ 
tion is significant. Hence, at present when no signatures 
of these new particles have been found and, as a result, 
the direct search experiments are regularly pushing the 
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lower limit of the allowed masses of these new particles 
upwards, the question gets increasingly important. 

In this paper, we want to illustrate that assuming 
single productions to be negligible compared to pair 
production may be unnecessary (or even improper) in 
some cases. In fact, including single productions in the 
theory estimation might give us extra clues about the 
BSM models being probed. We shall use the example 
of the recent leptoquark (LQ) search by CMS [1] to 
indicate how single productions can be systematically 
included in this search for pair production of scalar LQs 
and demonstrate how the extracted exclusion limits 
(ELs) on parameters (like masses, branching fractions 
(BEs) etc.) change in presence of model dependent 
production processes with similar final states as the pair 
production which could potentially ‘contaminate’ the 
signal simulations. 

LQs are hypothetical color-triplet bosons (scalars or 
vectors) that also carry lepton quantum numbers. Hence, 
a heavy LQ can decay to a lepton and a quark. They ap¬ 
pear in different BSM theories like the Pati-Salam models 
[2], SU(5) grand unified theories [3], models with quark 
lepton compositeness (see e.g. [4]), colored Zee-Babu 
model [5] etc. In some supersymmetric models squarks 
can also couple to a quark and a lepton via some R-parity 
violating couplings [6] and thus have similar phenomenol¬ 
ogy to the LQs. The leptons in their decay states make 
their collider signature relatively clean compared to that 
of other colored particles that dominantly decay hadron- 
ically. At the LHC, searches for LQs have been going on 
for quite some time. Both ATLAS and CMS have put 
limits on LQ masses [1, 7-13]. 

Recently, CMS collaboration has reported the results 
of their searches for the first generation scalar LQs at 
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FIG. 1: Representative parton level Feynman diagrams for pair [(a)-(e)] and single productions [(f)-(j)] of LQs at the LHC. 


the 8 TeV LHC with 19.6 fb“^ of integrated luminosity 
[1]. They have searched for pair production of LQs in 
two different channels - i) the eejj channel where both 
LQs decay to electrons and jets and ii) the eujj channel 
where only one LQ decays to an electron and a jet while 
the other one decays to a neutrino and a jet. They have 
found a mild 2A(j {2.6a) excess of events compared to 
the Standard Model (SM) background in the eejj {evjj) 
channel for LQs with mass around 650 GeV. The 

excesses observed certainly make the LQ search more in¬ 
teresting to be investigated further [14-19] . The airalysis 
also puts 95% CL ELs on first generatiorr scalar LQs for 
< 1005 (845) GeV for /3 = 1 (0.5), where (3 is the 
BF for a LQ to decay to an electron-quark pair. 

Here, we take a closer look at the CMS analysis and in¬ 
vestigate LQ single production processes that could give 
potentially significant contribution to the signal [20] but 
were ignored in the analysis. In the next section, we 
specify the single production processes under consider¬ 
ation, in Sec. HI we present the phenomenological LQ 
models and discuss the method that we use to compute 
single productions, in Sec. IV we present our results and 
finally, in the light of the new results we further clarify 
our argument and conclude in Sec. V. 


II. LQ SINGLE PRODUCTIONS AT THE LHC 

The limits from HERA data [21, 22] roughly indicate 
that the generic (scalar) LQ-lepton-quark coupling, A 
could even be larger than 0.5 if > 600 GeV. How¬ 
ever, the CMS analysis assnmes that A is small enough to 
ignore all LQ single productions and sets A = Aqed = 0.3 
to compute the total decay width of a scalar LQ for the 
signal simulations. In the leading order (LO) of cou¬ 
plings, the pair production of LQs is almost independent 
of A (see Figs, la - Id) except for the t-channel lepton 
exchange contribution whose amplitude is proportional 
to A^ (Fig. le). For small A, this contribution is actually 
quite small, e.g. with A = 0.3 it amounts to only about 5 
percent of the A independent contribution to the LO pair 


production cross-section for Mi^ = 650 GeV (see Table 
I). This is why it is justified when the A^ contribution is 
ignored in Refs. [23, 24] to compute the next-to-leading 
order (NLO) cross-sections for the pair production (used 
in the CMS analysis). However, while it may be safe to 
ignore the t-channel lepton exchange contribution to the 
pair production for A = 0.3, we shall see that even with 
this value of A, single production of LQs could actually 
contribute quite significantly to the analysis. 

The important point here is that a fraction of the sin¬ 
gle production events would pass the signal selection cuts 
used in the CMS analysis and contribute to their esti¬ 
mations of pair production events. For example, in the 
eejj channel, the analysis considers pp {£q iq) —>■ ejej 

(the curved connections above or below mark a pair com¬ 
ing from the decay of a leptoquark, iq) and hence de¬ 
mands the presence of exactly two electrons and at least 
two jets in the selection. However, for non-zero A, any 
pp —>■ {iqej) —>■ ejej process could also contribute to the 
signal. Hence, a priori, these should also be included in 
the signal simulations. Since the CMS analysis requires 
at least two jets for event selections, the following two 
types of hard processes with 3-body final state can con¬ 
tribute: 


1. Born single production with radiation (BRl): The 
Born diagrams for pp —>• iqi (as shown in Figs. If & 
Ig) can not directly contribute to the CMS search 
but with the emission of a QCD radiation (we refer 
to both parton splitting and radiation simply as 
‘radiation’) they contribute at O 


2. New subprocess of three-body single production 
(NS3): It originates at 0{g^X) and includes dia¬ 
grams for pp —^ iqij that do not count as BRl 
single production (see Figs. Ih-lj for example). 
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III. LQ MODEL AND SIGNAL SIMULATIONS 


As already stated, CMS assumes that a generic first 
generation scalar LQ decays to an electron and a quark 
with a BF /3 and a neutrino and a quark rest of the time. 
Hence, we also adopt a parametrization in which a LQ 
decays to a eq pair via a coupling Ag with = 

Pe = P and to a vq pair via with = Pv = 

(1 — P), and take 

xl = p,\^ = p\\ \i = p,\^ = (i-p)W ( 1 ) 


so that for any /3, the total decay width remains fixed at. 


i.e., the expression used in the CMS analysis. Motivated 
by the experiment, we consider scalar LQs that couple 
with only first generation quarks and take two simple 
Lagrangians for scalar LQs with electromagnetic charge, 
Qem = —1/3 (Model A) and 2/3 (Model B) respectively, 


■^int = (v^MRSl + 4 + Xyd^iVllq 

+ H.C., (3) 

~ (■\/A “t“ XyUYll'i^Eq 

+ H.C., (4) 


where rji^ and rjji = (1 — Tyn) are the electron chirality 
fractions, i.e., ryn (? 7 r) gives the fraction of electrons com¬ 
ing from a LQ decay that are left-handed (right-handed). 
In other words, r]i^Pe (??R/3e) is the BF for a LQ to decay 
into a left-handed (right-handed) electron and a quark. 
As the experiment is insensitive to the polarization of 
electrons, for simplicity and without loss of any general¬ 
ity, we have set rjL = 1 for our computation. 

The above models are completely generic and as long 
as no distinction is made between e'^ and e~, they can 
act as templates for a wider variety of LQs that can de¬ 
cay to SM quarks and leptons. For example, an analysis 


done in Model A for the eejj channel separately will also 
be applicable for a LQ with charge 5/3 that couples to 
a u-type quark and an electron (but not to a neutrino). 
Similarly, the LQ from Model B can represent a charge 
—4/3 particle also. Moreover, these models could be con¬ 
nected to the ones generally found in literature (see e.g., 
[25, 26]) via parameter rescaling. 

With these models, we have performed Monte Carlo 
(MC) simulations for the LQ single productions in both 
eejj and evjj channels (as long as the LQ width re¬ 
mains small, we ignore the interference between single 
and pair production). For our simulations we have em¬ 
ployed MadGraph5 [27] (and Pythia6 [28] within it) to 
generate and shower events and Delphes 3 [29] to sim¬ 
ulate the CMS detector environment and implement the 
selection cuts. We have used CTEQ6L1 Parton Distri¬ 
bution Functions (PDFs) [30] for all our numerical com¬ 
putations. 

A. The eejj channel 

Following the CMS analysis, we set P = 1 to gener¬ 
ate the signal events from the single productions in the 
eejj channel. Unlike the pair production, events from the 
lowest order single production process (strictly Born dia¬ 
grams) will not pass the hnal event selections as they will 
have only one jet in the final state. At the parton level, 
one could roughly estimate the (potential) contribution 
to the eejj channel coming from the inclusive single pro¬ 
duction (cr|®) by computing the cross-section for 

pp^ (iqej) -^ejej (5) 

without including the pair production diagrams. This 
will include contributions from both BRl and NS3 sin¬ 
gle production processes. But the BRl process contains 
divergent pieces that get cancelled when loop diagrams 
(virtual correction to Born diagrams) are properly in¬ 
cluded. The picture becomes clear if we express the 
cross-section for the inclusive single production with two 
electrons in the final state as. 


Divergent terms 


_ee _ f y-r^O 


-Virtual 


soft+collinear 

DRl(^qej') 


+cr- 


soft 

NS3(^gei) 


) + ( 


,hard 


^hard 


I IlCllU \ 

+ ^NSSie^ej)) 


( 6 ) 


Small contribution 


Main contribution 


where the neglected terms are of higher order than 
O {alX^) and the superscripts ‘soft -I- collinear’, ‘soft’ 
or ‘hard’ refer to the nature of the extra jet. Notice 
that, just like the ‘pure’ Born process, the pieces with 
divergence ultimately contribute very little to the exper¬ 
iment (if only parton showers are included) because it 
demands two separated ‘hard’ jets (see the selection cri¬ 
teria in [1]). Hence, in a tree level calculation, one could 


put some minimum p^-cut on the second jet (the jet com¬ 
ing from the LQ is generally the hardest one) to avoid the 
small-contributing part. This way one could compute the 
cross-section, but it would then become a function of the 
PT-cut and hence ambiguous (since a priori there is no 
absolute way to choose the ‘ideal’ p^-cut) and unsuitable 
for our purpose. 
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TABLE I: Cross-sections (in fb) for LQ productions at the 8 TeV LHC. The cross-sections are obtained with /tr = /rp = Mi^. 
The inclusive single production cross-sections, = /3e x (j(pp — >■ £qe + nj) « (see Eq. 8) and ctJ‘^(/3) = x cr(pp — >■ 

IqC + nj) + Pe x o-{pp IqV + uj) « 2A^/?(1 — P)^^'' (see Eq. 10) with n > 0 are obtained by ME©PS technique (as described 
in Sec III A) with Model A (Eq. 3). 


LQ 

Mass 

(GeV) 

Pair production 

0-p°(pP -f £qtq) 

A 0 A = 0.3 A = 0.5 

Inclusive single production 

gT{P = 1.0) 

A = 0.3 A = 0.5 

oTiP = 0.5) 

A = 0.3 A = 0.5 

550 

650 

750 

24.4 22.4 21.8 
7.2 6.6 6.4 
2.4 2.1 2.1 

41.2 118.7 

17.7 51.3 

8.3 24.5 

17.0 48.2 
7.3 20.5 
3.5 9.7 


TABLE II: MC events for LQ signals obtained for the 8 TeV LHC and 19.6 fb“^ of integrated luminosity. For the pair 
production, we quote the number of MC signal events obtained by CMS, J\fp [1] (with (t^'^°[ 23]). For single productions, we 
show the number of events obtained with Model A (Eq. 3) that survive the same final eejj or eujj selection criteria as the 
pair production for each (described in the CMS report). A/'s. We also show the selection-cut efficiency, {Mi^) - the 

fraction of single production events that survives the selection criteria for a particular Mi,^ (see Eq. 12). 


LQ 

Mass 

(GeV) 

eejj channel (A = 0.3, P = 1.0) 

eujj channel (A = 0.3, /3 = 0.5) 

Pair prod. 
Events (Afp) 

Single prod. 
Events (A4"Q el" (%) 

Pair prod. 
Events (A/j,) 

Single prod. 
Events (A4"") (%) 

550 

650 

750 

410.5±1.9 

125.9±0.6 

43.1T0.2 

83.0 10.3 

30.8 8.9 

12.4 7.6 

121.4±1.2 

37.2±0.4 

12.9T0.1 

18.1 5.4 
6.2 4.3 
2.6 3.8 


Clearly, a better approach would be to consider loop 
diagrams consistently to compute the pp —f iqe processes 
upto NLO QCD level (keeping the hard radiation) includ¬ 
ing the NS3 contribution to compute the inclusive single 
production cross-section (fT|®) correctly upto O (ctgA^) 
without any ambiguity. We are at present working in 
this direction [31]^ but for this paper we employ matrix 
element-parton shower matching (ME0PS) technique to 
obtain a theoretical estimate of the inclusive single pro¬ 
duction cross-section without any arbitrary pp cut on 
the second jet. The main difference between a proper 
NLO QCD computation and our ME©PS computation 
will come from the part that contributes very little (from 
the second and third terms of Eq. 6). Hence, for the 
same luminosity, the number of events passing the se¬ 
lection criteria (require atleast two jets with pp > 45 
GeV) [1] should not vary too much between a proper 
loop included NLO level computation and the one with 
ME©PS. The estimates for the inclusive cross-section will 
be different but this difference should not matter if the 
efficiencies of the selection cuts (will be defined in Sec. 
IV) are computed consistently. 


Eor the ME©PS computation we have used the shower- 
kp scheme [34] . We generate events for the inclusive sin¬ 
gle iq production signal as the combination of the follow- 


^ Although not directly useful for this paper, we note that some 
estimates of the NLO LQ single productions at the 14 TeV LHC 
in different models already exist in the literature (see e.g. [32, 
33]). 


ing processes. 


pp 

pp 

pp 


{£q e) e , 

ej) -f ej , > 

i^q ejj) -f ej ejj . ^ 


( 7 ) 


The A and /3 dependence of inclusive single production 
cross-section in the eejj channel can be made explicit as: 

ar(/3,A,M,J = A2^eSr(M,J+... 

(MO+---, (8) 

where the neglected terms are at least 0{X'^). Notice 
that with the chosen parametrization, the P dependence 
of the single production cross-section becomes the same 
as the pair production. 


B. The evjj channel 

For the evjj channel we set ,0 = 0.5 to generate the 
signal events with a similar ME©PS computation (the 
argument for computing this manner remains same as 
before). 


pp 


i^q 

e 

-b 

^q 



ejj^, 


pp 

-)► 

(^q 

ej 

-b 

iq 

vj) 

-f 

ejjp j, 

> (9) 

pp 

-)► 

(^q 

ejj 

-b 

iq 

vjj) 

-)■ 

ejjp jj ■ 



Like the eejj cross-section, here too the A and (3 depen¬ 
dence of the inclusive single production cross-section can 
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Mi, (GeV) 

(a) eejj channel (taken from [1]) 



Mt, (GeV) 

(c) eejj channel (Model A: QemC^?) = —1/3, 5/3) 



Me, (GeV) 


(b) eujj channel (taken from [1]) 



Me, (GeV) 

(d) evjj channel (Model A: Qbm)^?) = —1/3)) 


FIG. 2: Effect of single productions on exclusion limits. On the left panel: eejj channel and on the right panel: evjj channel. 
See Eqs. 8 & 10 for the definitions of dl® and respectively. The red solid lines are theoretical estimations. The upper-row 
plots (A —>■ 0) are same as the ones presented in Ref. [1]. Eor the lower-row plots (A = 0.3, the value used in the CMS analyis), 
the expected and observed upper limits at 95% CL on the LQ production cross-sections are obtained after correcting by 
(A, J (scc Eqs. 15 & 19). 


be made explicit as, 

ar 

= \l /?. {M(,) + Xl He (M,J + • ■ • 

= {Ml,) + A^ {Mi,)+--- 

A^ 2^(1-^) dr , (10) 

where, like before, the neglected terms are at least 0{\^) 
and 

dr {Ml,) = i (dr {Ml,) +dr (m,j) . (n) 

Again we see that with the chosen parametrization, the 
P dependence of the single production cross-section be¬ 
comes the same as the pair production. 


IV. RESULTS: RESCALED EXCLUSION LIMITS 


In Table I we show the theoretical cross-sections for 
pair production at the LO and inclusive single produc¬ 
tions (in Model A) for eejj and evjj channels with Mi = 
550, 650, 750 GeV and A = 0.3, 0.5. As expected, the 
pair production is almost independent of A whereas the 
single production cross-section increases as A^. In Table 
II we show the number of events that pass the final se¬ 
lection criteria for pair production (A/p) (from Ref. [1], 
obtained with o-p [23] ) and inclusive single production 
(A4) for both eejj or evjj channels. In the same table 
we also show the selection-cut efficiency for the single 
production, Cg - the fraction of events that survives the 
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Jl, 



M,, (GeV) 

(a) Model A: QEM(^g) = —1/3, 5/3 



Jl_ 



Mt, (GeV) 

(c) Model B: Qbm(^(j) = 2/3, —4/3 




FIG. 3: The exclusion limits (solid or dashed lines) from the observed data in A — Me^ plane for // = 1.0 (left pannel) and 
j3 = 0.5 (right pannel) for Models A and B. All points on the left of the solid or dashed lines (shaded regions) are excluded by 
the data at 95% CL. When f) = 0.5, the strongest limits come from the evjj data for both models except in Model A the eejj 
data gives the strongest limit for A > 0.55 (M^^ > 910 GeV). 


selection criteria for a particular , 

where C = 19.6 fb“^ is the integrated luminosity for 
which Afs is computed and X = {ee or ez/}. Note that 
does not depend on any overall factors in Us like 
A^,,5g (for eejj channel) or jSePv (for evjj channel) etc. 
The numbers show that even for A as small as 0.3 a good 
number of events survive the selection cuts optimized for 
the pair production. Moreover the ratio of Afs/N'p in¬ 
creases with increasing i.e., with increasing mass, 

single productions contribute more and more compared 
to the pair production. 

In Fig. 2, we show how the presence of single produc¬ 


tions affects mass ELs. In the vertical axes of these plots 
we have now plotted 

/3^(Tp-I-(^, A) =/3^ (cTp-h A^dg®) (13) 

instead of /3^crp for the eejj channel and 

2/3 (1 - /3) ap + (/3, A) =2/3(1- /3) (ap + 

instead of 2/3 (1 — /3) CTp for the evjj channel where 
and al" are defined in Eqs. 8 & 10 respectively. When 
A —> 0, single productions vanish, hence the upper row 
plots are identical to the ones presented by CMS (we 
avoid writing A = 0 as in that case both the decay chan¬ 
nels would vanish). For the plots with A = 0.3, we have 
rescaled the expected and observed 95% CL upper lim¬ 
its on the LQ pair production cross-sections obtained by 
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1 

0.8 

0.6 

0.4 

0.2 



300 400 500 600 700 800 900 100011001200 


A4 (GeV) 

(a) eejj channel (Model A: Qbm(^ij) = —1/3, 5/3) 



Mt, (GeV) 

(b) eiyjj channel (Model A: QEM(^q) = ~l/3) 



Mi^ (GeV) (GeV) 

(c) eejj channel (Model B: Qem)^?) = 2/3, —4/3) (d) eujj channel (Model B: Qem(^ij) = 2/3) 

FIG. 4: Regions excluded at 95% CL by the observed data in the P — X plane for different values of A in Models A and B. The 
plots for A —1- 0 are obtained from pair production only and hence are same as the ones presented by CMS [1]. 


CMS by multiplying them with a factor 


and Cp+s is defined as follows, 


TZ^{X,MeJ = 




■"p+s 




(15) 


where Cp is the selection-cut efficiency for pair production 
defined in the same manner as eg in Eq. 12, 


= 


(/X(^)^NLO) X/: ’ 


with 


(/3) = for X = ee , 

= 2/3(1 — /3) for X = ev 


(16) 


(17) 


Ep+s 

(/3) cTp (M,J (M,J + gX ef (M,J 


/X (/3)ap(M,J+aX(/3,A,M,J 
^ gp (M,J (M,J + A^df (M,J ef (M,J 

CTp (MfJ + A^ffX 

When A —0, [P, A, vanishes and hence 

lim TZ^ iX,Mi ) = 1 . 

A-s-O ^ ' 


(18) 


(19) 


For our analysis, we have used the numbers given in the 
CMS analysis to compute ep. 
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(c) eejj channel (Model B: Qem(^(i) = 2/3, —4/3) 



(d) eujj channel (Model B: Qem)^?) = 2/3) 


FIG. 5: Regions excluded at 95% CL by the observed data in the /3 — A plane for different values of in Models A and B. 
On the left panel: eejj channel and on the right panel: evjj channel. The exclusion limits for fixed Mt^ values are shown 
as boundaries between differently shaded regions (except the Miq = 650 GeV lines which are drawn explicitly). For the eejj 
channel, the region above any fixed mass line is excluded whereas for evjj channel the region on the right of any fixed mass 
line is excluded. 


We can see how the 95% CL ELs change with A from 
Fig. 3 which shows the observed ELs in the A — 
plane for /3 = 1.0 and fj = 0.5 and both Models A and B. 
These plots can also be used to put upper limits on A (or 
equivalently on Ag and Ai, for any /3) for certain values of 
Mi^. Once single productions are considered, the limit 
coming from the eejj data for A = 0.3 (0.6) on the mass 
of the LQs with Qem = —1/3,5/3 goes up to about 1070 
(1230) GeV from 1005 GeV for fj = 1.0. However, for 
the LQs with Qem = 2/3, —4/3, the corresponding limit 
reaches only about 1030 (1110) GeV. When /3 = 0.5, the 
most stringent limits come from the evjj data for both 
models except in Model A where the eejj data gives the 


strongest limit for A > 0.55 {M(^ > 910 GeV). For a LQ 
with Qem = —1/3, the mass EL improves from 845 GeV 
to about 870 (970) GeV for A = 0.3 (0.6). For the LQs 
with Qem = 2/3 the same limit improves to about 880 
(930) GeV. The excluded regions in the f3 — Me^ plane are 
shown in Fig. 4 for A ~ 0, 0.3 and 0.5 respectively. The 
plots show that for any fixed /3, the ELs increase with 
increasing A. In Fig. 5, we show how the observed ELs 
vary with /3 in the /3 — A plane for the two channels in 
both Models A and B. Notice, that in the evjj channel 
plots (Figs. 5b & 5d), the EL curves for = 650 GeV 
lie among the = 700 and 800 GeV curves. This 
happens because of the observed excess at 650 GeV in 
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700 800 900 1000 1100 1200 1300 1400 


M,, (GeV) 

FIG. 6: Comparison of 95% ELs obtained by including the 
contribution of single productions in the pair production 
search data [1] (single in pair) and that of the pair production 
in the single production search data [13] (pair in single). Here 
we have considered Model A to simulate the single production 
events and set /I = /3e = 1. 


the data. 

Inclusion of single productions in the analysis can 
somewhat improve the x^-fit of the observed data in eejj 
and evjj channels separately. We use the following 
function, 


(A A) 




where IVlq, .A^bg and IVoata are the number of LQ sig¬ 
nal events (MC), SM background events [1] (MC) and 
observed events (data), respectively for the bench¬ 
mark LQ mass taken by CMS (varied from 300 GeV to 
1200 GeV in steps of 50 GeV). The errors, AiVgQ are 
obtained by adding the systematic and statistical uncer¬ 
tainties of the total background in quadrature. In Model 
A with = 650 GeV and A = 0.6, as we move from 
‘only pair’ to ‘pair-hsingle’, the best fitted /? changes from 
0.19 (0.044) to 0.15 (0.025) in the the eejj {evjj) chan¬ 
nel. The corresponding decrease in the minimum of 
{Axmin) about 18% (15%). Similarly, for Model B, the 
best fitted /3 changes from 0.19 (0.044) to 0.17 (0.022) in 
the the eejj {evjj) channel and the corresponding 
is about 13% (11%). However, when we combine the two 
channels Ax^in becomes insignificant in both the mod¬ 
els. 

Before we move on, we make a small digression here. 
Very recently, GMS has presented a new dedicated search 
for single production of scalar LQs with a final state of 
two electrons and one jet (for the first generation) [13]. 


If we compare roughly, for Model A with A = 0.6 and 
/3 = 1, our rescaled EL from the pair production data on 
(1230 GeV) is competitive to the one obtained (1260 
GeV) there. For A = 0.4 and (3=1, the rescaled limit 
actually gives a better estimation of the lower bound on 
than the single production, because the former one 
includes the contribution from pair production. Now, 
just as we have included the contribution of single pro¬ 
ductions in the pair production search, we can apply the 
same logic and consider the contamination from the pair 
production in this dedicated search for single production 
by similar rescaling. In fact, in this case one can not 
avoid the pair production completely by tuning A, since 
it is almost independent of this parameter. In Fig. 6 we 
compare the limits thus obtained from these two cases 
in the A — plane. The EL obtained by including 
the single production contribution in the pair production 
search (i.e. same as in Fig. 3a) is shown by the solid 
line and the one obtained by doing vice versa is shown 
by the dashed line. For A < 0.55, the pair production 
search data gives the better limit on M(^. However, with 
increasing A and (excluded) single production con¬ 
tribution increases and pair production contribution goes 
down. For A > 0.55 and > 1200 GeV, the single pro¬ 
ductions take over the pair production to determine the 
EL. For such large A and , the contribution of the pair 
production becomes negligible and hence, in this range, 
both limits converge to the one obtained by considering 
single productions only [13]. 


V. DISCUSSIONS & CONCLUSIONS 


One must be careful while interpreting our results since 
they are obtained by simple rescaling, not a full statisti¬ 
cal analysis. Nonetheless, they certainly show that within 
the probed range of , systematic inclusion of the sin¬ 
gle productions in the signal simulation of the pair pro¬ 
duction search improves the ELs. For A = Aqed = 0.3 
(the value used for MC simulations) the CMS analysis 
in Ref. [1] underestimates the mass ELs. This is just a 
manifestation of the fact that the reach of the LHC (to 
probe new particles) can increase if multiple production 
channels are systematically combined (we have already 
demonstrated this while estimating the LHC discovery 
reach for leptogluons which also decay to £j pairs like 
LQs [35], support of this argument can also be found in 
Ref. [20]). While it can be argued that the pair pro¬ 
duction gives the most conservative estimation of ELs on 
mass of LQs (or any colored particles), but, as we have 
seen, it can be misleading about other parameters like 
the branching fraction, (3. Moreover, by ignoring single 
productions completely, one also ignores the fact that 
the pair production data is not only excluding masses of 
new particles but also putting limits on new couplings. 
For LQs, the limits on A thus obtained are competitive 
with the ones obtained from HERA observations [21, 22]. 
Our naive estimates indicate that compared to the limits 
from HERA, the CMS data already provides better lim- 
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its (both in terms of lower limits on and upper limits 
on A) for the entire parameter range covered in Fig. 3. 

The results also support our argument that for any 
fixed A, the ELs depart more and more from the pair 
production ones as increases. Hence, if single pro¬ 
ductions must be ignored, it also has to be made certain 
that within the whole mass range probed, the single pro¬ 
ductions remain small enough and the choice of A should 
reflect this. Alternatively, the selection-cuts could be op¬ 
timized in such a way that they disfavor single produc¬ 
tions strongly. But, rather than making any assumptions 
like this, our suggestion is to include single productions 
in the signal simulations and allow A to vary as much as 
the existing experimental limits allow. Moreover, if the 
selection cuts could also be optimized so that they do 
not disfavor the single productions unnecessarily, then 
it may be possible to impinge more on the presently al¬ 
lowed parameter space due to increased sensitivity to A 
(also see e.g. Refs. [20, 35]). Also, advanced techniques 
(like multivariate analysis) could tell us more about pa¬ 
rameters other than the mass. This way, from the same 
experiment, we might learn more about the underlying 
model than a conservative mass exclusion of LQs. 

So far our discussions have been mainly centred on the 
example of LQs. However, it is easy to see that the ar¬ 
gument in favour of including single productions is appli¬ 
cable for many other BSM particle searches. Generally, 
if there is no symmetry that prevents model dependent 
production (like R-parity prevents single production of 
supersymmetric particles), ideally all such productions 
that could (potentially) contribute to their search should 
also be considered. At the core of our argument stands 
the very simple and well known statement - ideally^ all 
processes that can have similar final states should be in¬ 
cluded in the signal simulations. However, in practice, 
one often has to make simplifying assumptions like, if 
there are too many processes to produce a new parti¬ 
cle that could contribute to its search, usually the ‘sub¬ 
dominant’ ones are ignored, or if there are different new 
particles that can have similar signals, generally only one 
is considered at a time etc. Here, we are simply trying to 
suggest that when making such simplifying assumptions 
one should justify them quantitatively to the extent pos¬ 
sible and, moreover, if technology permits one need not 
make any ‘unnecessary’ assumptions. For example, in 
the context of LQs, for very large (say ~ 1500 GeV) 
there is no experimental bound that prevents A from be¬ 
ing order one. Hence, while probing such masses one need 
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not set A to be very small. Similarly, as our estimation 
suggests, a dedicated search for single production of LQs 
at the LHC should also consider the possible contribution 
from pair production. 

To summarize, we have shown that while searching for 
new colored particles at the LHC, one need not always 
consider model dependent (single) and (mostly) model 
independent (pair) productions separately. It may be un¬ 
necessary to ignore single productions compared to the 
pair production in some cases. Single productions can 
obviously be important if the couplings controlling them 
are not small - this is why the studies that focus on 
single productions generally consider large couplings (see 
[36] for an example in the context of LQ) - but even if one 
assumes the couplings to be small, single productions can 
still contribute significantly if the mass of the particle be¬ 
ing probed is sufficiently high. In such a case, by ignoring 
single productions naively, one might conclude wrongly 
about parameters like branching fraction etc. With the 
example of the recent search for pair production of scalar 
LQs performed by GMS, we have illustrated how the con¬ 
tribution from single productions can be included in the 
analysis and how it can lead to modified limits on the LQ- 
mass. In addition to this, inclusion of model dependent 
single productions gives us new information about model 
couplings like the LQ-lepton-quarks couplings from the 
same experimental data. For this purpose, we have pro¬ 
posed two generic models for scalar LQs that can act as 
templates for a wider variety of LQs and can accommo¬ 
date the experimental data easily. Applying similar logic, 
we have also demonstrated how the pair production could 
affect bounds obtained from the more recent GMS search 
for LQs in single production channels. Finally, we have 
pointed out that our arguments are not specific to the 
case of LQs (which we have used as an illustrative ex¬ 
ample), but they are applicable to other BSM searches 
too. 
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